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ABSTRACT 


A simple but computationally efficient finite difference 
scheme has been used to solve the problem of conjugate heat tran- 
sfer in a laminar flow of high Prandtl number fluid in a circular 
tube whose periphery is exposed to non-uniform radiation heat 
flux from a large heated wall. The present model is a simulation 
of the effect of hot refractory wall close to a tube in a pipe 
still. 


The results Indicate that in a thin tube with low conduc- 
tivity, the temperature of the tube is very high which may result 
in hot spots on the tube surface* If the tube material is thick 
and has a sufficiently high conductivity, the peripheral conduc- 
tion in the tube wall alleviates this problem. Also, the temper- 
ature of the tube goes up when the tube is placed closer to the 
heated wall. 

For the wider wall the mean temperature of the inner 
wall of the tube is more as the tube receives more radiation* 
Also, it is observed that mean tube inner wall temperature incr- 
eases with the axial distance and is lower for higher Prandtl 
number fluid. 

It is also seen that thermal entry lengths are higher 
for higher Prandtl number fluid and for higher Reynolds number 
of the flow. 

The fluid bulk temperature increases with the axial 
distance and for higher Prandtl number of the fluid and lower 
Reynolds number of the flow, the fluid bulk temperatures are 
less . 


The accuracy of the present numerical solution has been 
checked by comparing it with the analytical solution of the case 
of fully developed constant heat flux condition on the periphery 
of a very thin tube* The result is found to be very satisfactory 


CHAPTER 1 


INTRODUCTION 


General 

This investigation presents the case of heating of a 
fluid flowing through a long circular tube which is exposed tc 
a non-uniform radiation heat flux on its periphery by the pro) 
mity of a very large heated wall situated parallely above the 
axis of the tube. The temperature of the fluid is unifoim at 
the inlet of the tube where the heat transfer begins, but the 
velocity profile is already fully developed and invariant. 
This, however, limits the applicability of the present analyst 
to fluids whose Prandtl numbers are very high relative to 1 , 
for example, oils. The present analysis is for constant 
property fluids but variable property fluids can also be inve< 
tigated with little modification in the solution procedure. 

The present model is a simulation of the effect of ho 
refractory wall close to a tube in a pipe still. Frequently, 
in the case of uneven distribution of heat flux around the 
tube, hot spots on the tube surface may occur. If the tube 
material is thick and has a sufficiently high conductivity, 
peripheral conduction in the wall may alleviate . this problem, 
but for thin walled tubes, the difficulty can be accute. 

The temperature distribution in the wall is a strong 
function of the amount of radiation heat flux received at the 
outer surface of the tube and the heat transfer at the interf 
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between the inner surface of the tube and the fluid. Thus, 
the conjugate heat transfer at the interface plays an impor- 
tant role in this problem, 

1.2 Objective 

The objectives of the present investigation are (a) to 
find the temperature distribution in the wall of the tube for 
various thickness and thermal conductivity of the tube and for 
various positions of the tube with respect to the heated wall, 
(b) to get a qualitative indication of the thermal entry length 
for fluidsof different Prandtl numbers and for various Reynolds 
number of the flow, and (c) to show circumferential variation 
of solid-fluid interface temperatures for different widths of 
the heated wall, 

1 . 3 Literature Survey 

The problem of laminar flow in a circular tube with 
fully developed constant heat rate conditions axially, but with 
an arbitrary peripheral variation of heat flux, has been treate 
by Reynolds [1, 2], In [l ], a solution is obtained for a pulse 
of heat transferred across an elemental section of wall, and 
superposition is employed to obtain solutions for any prescribe 
heat flux. In [2], the prescribed heat flux is expressed in 
terms of a Fourier expansion and also the flow is turbulent. 

An interesting conclusion of [2] is that the effects of circum- 
ferential heat flux variation in turbulent flow are sometimes 
more pronounced than in laminar flow. The coupled effect of 
wall and fluid conduction in laminar pipe flow has been studied 
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by Faghri and Sparrow [ 3 ] and Zariffch et al [4] using finite 
difference and by Gampo and Rangel [b] analytically. Davis 
and Venkatesh [6] also studied heat transfer to Poiseuille 
flow in a pipe with heat conducting section, as an example of 
their new general formulation of the solution of conjugate 
boundary value pi^blem. Computations, however, were not per- 
formed. Luikov et al [ 7 ] presented a closed form analytical 
solution for the conjugate heat transfer problem in case of 
fluid flow in circular tube* S.V. Patankar [8] suggested the 
use of harmonic mean of the transport properties on both side 
of an interface to treat the points of discontinuity. Therm- 
ally fully developed flow in a square duct with finite wall 
thickness was presented as an example. The method is quite 
general and can also calculate fluid flow and heat transfer i 
arbitrary geometries. Barozzi and Pagiliarini [9] have used 
finite element method to solve axi- symmetric conjugate heat 
transfer in a pipe with uniform heat flux specified on its 
periphery. The method combines a finite element treatment ol 
the energy equation in the tube wall and the application of 
Duhamel' s theorem at the interface. 

1 *4 Salient Features of the Present Investigation 

The present investigation differs from the previous 

works on the following counts* 

(i) the non-uniform peripheral heat flux on the outer 

surface of the tube is actually calculated, instead 
of being prescribed. 



(ii) 
(iii ) 
(iv) 


(v) 


the flow is thermally undeveloped- 
the problem is non-axi symmetric, 
a simple but computationally efficient finite 
difference scheme has been used to solve the 
problem. 


it has the capability of handling variable 
property fluid and tube material with little 
modification in the numerical method. 
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CHAPTER 2 

PROBLEM FORMULATION 
2.1 Description of the Problem 

The physical model is depicted in Figure 1(a) and 
Figure l(b). The refractory wall has a width W. The pipe, 
which has inside and outside radii r^ and r2, is separated by 
a distance d vertically from the wall, the tube being at a 
distance L from the left extreme point of the wall. L and d 
together determine the position of the tube with respect to 
the wall. The wall behaves as a black body radiator at tempei 
ature The surrounding is at temperature Tg, The upper 

portion of the periphery of the tube receives radiation heat 
flux both from the wall and the surroundings while the bottom 
portion receives radiation heat flux primarily from the surr- 
oundings. So, there is non-uniform radiation heat flux on the 
periphery of the tube (Figure 1(a)), but the radiation heat f] 
is uniform in the axial direction since both the wall and the 
tube are very long in that direction (Figure 1(b)). A fluid 
enters the tube at a uniform inlet temperature T^ while the 
velocity profile is fully developed. 

So, the thermal boundary layer develops along the 
axis of the tube until the temperature profile is fully deve- 
loped. The local heat flux density at any x position from th« 
inner pipe wall to the fluid is, 

q, = q,(0) = h(T, - Tjj,) 

where T^j^ is the fluid bulk temperature at that x-position 
and h is the heat transfer coefficient at the given x (the 



w 



FIG. 1 (a) SCHEMATIC DIAGRAM OF THE PROBLEM AND TH 
COORDINATE SYSTEM IN (r,c^) DIRECTION 



FIG. 1(b) SCHEMATIC DIAGRAM OF THE PROBLEM AND THE 
COORDINATE SYSTEM IN (r,x ) DIRECTION 
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method of determination of their values shown later). T. is 

1 

the wall temperature at same x- station. The polar coordinates 
(r, 0 ) are shown in Figure l(a). 

The corresponding radiant flux density to a point on 
the outer surface where the temperature T^ = 

q2 = qgCia) = - T^) 4 (1 - . T^)] 

The geometric view factor is the fraction of radi- 

ation leaving a small area on the outer pipe surface that is 
directly intercepted by the wall; and the remaining fraction, 

^ ^ 2W absorbed by the remainder of the surroundings. It is 

described in details in Appendix A, 

The temperature distribution in the wall of the tube 
and the fluid at a particular x-location are determined by sol- 
ving simultaneously the energy equation in the wall and the 
fluid and the interface and then the solution marches in the 
x-di recti on . 

Following assumptions are made to model the problem 
matheraati cal ly s 

Ca) The heated wall behaves as a black body radiator. 

(b ) The surrounding atmosphere is transparent to thermal 
radiati on . 

(c) Matural convection from the heated wall to the tube 
surface is neglected, radiation being the principal 
heat transfer mode. 

(d) The flow is steady, laminar and incompressible, 

(e) Physical properties of both the fluid and the tube 
wall are constant. 
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(f) The tube is smooth so that frictional loss is negligible. 

(g) Axial conduction in the fluid is neglected as in the 
present case RePr is greater than 100. 

(h) Axial conduction in the tube wall is neglected as the 
tube is sufficiently long. 

(i) Viscous dissipation is neglected as it is a low speed 
flow. 


2.2 Sasic Equations 

For incompressible, constant property fluid, the comb- 
ined momentum and continuity equation for fully developed flow 
becomes, 


1 d du\ 1 dp 

r dr dr ** p dx 


0 


with the boundary conditions, 

(i) no slip condition, u = 0 at r = 

(ii) the velocity is maximum at the centre i.e. ^ = 0 at r = 0. 
Solving the momentum equation subject to the above boundary 
conditions, the velocity profile obtained as follows! 


u(r) = - ^ if 

u = Average velocity (axial) 

_ %ax _ 1 ^ 

2 " §p dx 1 


as, 


u 

max 


u (r = 0) 


1 - ^ 

4p dx 1 


Inserting u in the expression of u(r), the velocity profile 
gets the form 


9 


( 2.1 ) 
( 2 . 2 ) 

0 (2.3) 

at r = fr ~ ^ " radiation heat 

flux falling on the upper surface 

. , _ _ _ q^{0) _ ^ “ ^2^ + n - F^^(^S)} 

(T^ - T^)] 

•■• 3tr = r 2 , li = [Fav Tw + (1 - F 2 w)t| - rf] (2.3a) 

Where = configuration factor, determined by Hottel- 

Gross- string method (see Appendix A) 

■^2 =Ur=r2 = ^^^> 

(ii) Fluid region; 

Energy equation; u |^ = <Xf[^ + r H ^ 

(2.4) 

with boundary conditions, at r = 0, |-^ = 0 

0 3 ? 

at X - 0, T = T 

^ e 



u(r) = 2u u^^(r) 


where 


u#(r) = t - (^) 

^1 


Energy equations! 


(i) For solid region; + 1 _ 


ar^ ^ %2 3^2 


with boundary conditions. 



At the interface, the ecfuality of temperature and heat flux 
exists, so that, 
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at r 

and 



s 9r 



T 


. dj 
■f 3r 


(2.5) 


2.3 Fluid Bulk Temperature 

For the calculation of the convective heat transfer 
coefficient h(x) from the wall to the fluid along the axis of 
the tube (x-direction ) , it is necessary to define a character- 
istic temperature of the fluid commonly known as ’bulk tempera- 
ture’ or 'mixing cup temperature*. The reason is that in a 
tube flow there is no easily discernible free-stream condition 
as is present in the flow over a flat plate* For most tube or 
channel flow heat transfer problems the topic of central 
interest is the total energy transferred to the fluid in eithe 
an elemental length of the tube or over the entire length of 
the channel. At any x position, the temperature that is indi- 
cative of the total energy of the flow is an integrated mass- 
energy average temperature over the entire flow area. The 
numerator of equation (2.6) below, represents the total energy 
flow through the tube and the denominator of the same equation 
represents the product of mass flow and specific heat integ- 
rated over the flow area. The bulk temperature Is thus repre- 
sentative of the total energy of the flow at a particular 
location - (x) . 

Mathematically the bulk temperature for the fluid is 


defined as, 



1 1 


fb 


Thus 


■fb 


2iz 1 

/ / PC u T r dr d0 

0 o 5 

2-!t ^1 

/ / P C u r dr d0 

o Q " 


Tf^(x) 


( 2 . 6 ) 


The heat flux from the wall of the tube to the fluid at a given 
x-location and circumferential location 0, 0) is defined 

as. 


qi{x, |3) = -[- |I 


5 = Xslij (2.7) 

r=r.|,s ‘r=r^,s 


and the average heat flux at a given x is defined as, 


2 % 


^ ‘’i 


(2.7a) 


The average wall temperature of the tube at a given x is expres- 
sed as, 


21C 


T,(x) = T, d0 


( 2 . 8 ) 


where T^ = T'.jCxr 0) - Interface temperature = {2.8a) 

Local heat transfer coefficient, h^ can be found by the foll- 
owing heat flux balance equation; 


or, 


q^(x) - hjj(T| - T^jj ) 

(x) 

h ^ ! 

^1 - ^fb 


(2.9) 


The local Nusselt number, Nu is f.ound by, 
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Nu^ 


h^C2r,) 


( 2 . 10 ) 


ft'* 

r * 
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CHAPTER 3 
METHOD OF SOLUTION 


3. 1 Finite Difference Approximations 

The energy equations of the solid and the fluid along 

with the compatibility conditions at the interface are discre- 

tised using finite difference approximations as described below, 

A typical orthogonal mesh system for using the Finite 

Difference approximation is shown in Figure 2. 

The following finite difference approximations are used. 

o 

3 3 ^ 

(i) Using central difference formulae for the terms — k, ^ 

2- Sr 

and at the point (i, j), 


ar 


3^T 

3^T 


3T 


CAr)2 


- 2^1, J - ViJ 

1 


- 2Ti,j - 

(A0)^ 

1 


- 

2(Ar~) 


JJr 


where at = uniform step size in r-direction 
A0 = uniform step size in ^-direction 
i s= index in r-direction 
j = index in 0-dlrection. 

(li) Using backward difference formula for the convection term 

” at a x-station, k+1, 

ax ’ 

ax Ax 
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where T. . . = Temperature of the point (i, J) at the previous 

^ » j » 

jc- station . 

3 T 

Also term can be discretised by backward difference formula 
as. 


3r 


i+1 , j 


T- , ^ — T, 

Ar 


3»1«1 Discretization of the Energy Equation for Solid 

The equation (2.3) subjected to the boundary condition 
(2.3a} can be discretized as follows, 

Equation (2.3); r ^ ^ 

3^2 r 3r ^2 3^2 


Boundary condition (2.3a); at r = r^; 


ar 


oe 

k. 


[F 


2W 'W 


4 - 


(1 - 


- 



Refer to Figure 2. 


Tot, = Number of grid points in r-direction for the entire 
fluid-solid region 

= Number of grid points in ^-direction, both for solid 


and fluid region 

N - = Number of grid points in r— direction, only in the 
rr 

fluid region. 

As x-direction is taken infinitely long, no specific number of 
grid points are taken for this direction. By marching procedure 
the solution can be continued upto any length and it can be 
terminated as per requirement. 

The uniform, step sizes are. 


A 
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(rj - r,) 


ii0 ^ 


2 ti 


H 


and Ar. 


0 




Discretizing equation (2.3) by finite difference approximation, 
as 


1 1 


r? (ii0) 


2 ^ ^^i,j 


] 0 


Simplifying the above equation, 


Ar 2 


Ar 2 , 

+ ( 2 r 2 + r. ^ 2(^) 

Writing the above equation in the standard tridiagonal matrix 


form, 




(3.1 ) 


Where 


''i 


B. 

r 


2^i - ^i 

o Ar 2 


^i 


2rf + r. Ar 
i ^ 1 s 

- IT, .3^1 + 
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The above equation is valid for any interior point (i, j ) of 
the tube wall except the boundary points, i.e, for all nodes 
numbering from i = +1 to i = - 1 and in 0 -direction 

all points except the points at 0 = 0 and 0 = ( 2 ti - 0 ).i.e. 

j = 1 and j = N^. 

More correctly saying, the equation (3.1) is valid also 
for j = 1 and j = but for computational purpose the equa- 
tion is slightly modified for all points with j = 1 and j = N^, 
as described below. 

For j - 1, the point (j-1 ) physically represents the 
point 'M 0 ’'* Therefore the equation (3.1) is modified by repla- 
cing (j -1 ) by for j = 1 and (j +1 ) by 1 for j = N 0 , both 
changes occur in the expression of D^* 

Rewriting equation (3.1) for different j's the general 
form becomes, 


A. I. . . + B.T. . + C.T. . , = D. 

i i-1,j 1 1,3 11+1, 3 1 


(3.1 ) 


with 


A. = 2rr - r. Ar 
i — i s 

Bi = - [4 r2 + 4(^)^ ] 


for all 3 


and 


G. = 2r7 + r . Ar^ 

1 1 1 s 

Ar 2 

D. = - 2(t#) j for 2 ^ 3 < (N ^- 1 ) 


‘A0 


j+1 


^0 


2 

= - for j = 1 

Ar 2 

= - 2(ir^ fh.l A.N^-1^ 3 =N0 



J^T.^ - 
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Ftowever the equation (3.1) is only valid for interior points 
of the tube wall. 

To get the discretised energy equation for the upper 
surface, i.e, i = (r = r^), the boundary condition (2.3a) is 
used to modify equation (3.1) as described below. 

Boundary condition (2.3a) is, at r = r^ (i.e. i = N ), 


where 


^ ^ Tf ^ (a f 1 1 

3r " k ‘-^2vV - ^2W-' “■ 2 

S’ 

= temperature of the tube at r = r^, this is repre- 
sented by T* 

Putting T in place of T 2 and discretising equation (2.3a) 

4 


2Ar 

■^s 


IJ H. (1 - 


but as . term introduces non-linearity in the discretized 

i I j 

equation it Is necessary to linearize this term using Patankar* s 
approach [ 10 ] to linearization of the source term. 

According to this approach, -the non-linear term should 
be linearized in such a manner so that fast convergence is 
attained. It also suggests that when the coefficient of the 
transformed linear term is equal to the slope of the original 
nonlinear term the convergence is fastest, however this is not 
true always. 

Using this approach if . term is linearized as, 

1 j J 


where 



= Temperature at point 



j) in the previous 


iteration. 
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Putting this in the discretized equation of the boundary condi- 
tion (2.3a), the equation becomes, 

— ^ fi: A, ( ^ p \t4 

2Arg - ^ + (1 - 




Simplifying it, 


T. 


2cre 


i+i,j ' ~Tl f'^avKj^w + (' - ^aw >’’s + 

- T._j] + 

Putting this expression in equation (3.1) the discretized energy 
equation for i = is obtained as follows* 


hh-ij + 


(3.1a) 


where A. 

1 


B. 

1 


where 


D. 


"^^i " "^^2 ^ ^r ^i = ^2^ 

r, O /J r~ 2 

- - 4-2- 


= 0 


2 




- SjEFjiy Tyj + (1 - F2 w3'^S * 

J J 

2oe 


®2 = -ff + -2 ^-5^ 


The above equation is valid for i = only, here also 
is modified for j = 1 and j = as before, i.e, replacing 
(j-1 ) by for j - 1 and (j+1) by 1 for j = 
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3.1*2 Discretization of the Energy Equation for Fluid 

The energy equation for fluid (Equation (2.4)) is 
changed by inserting the expression of u(r) into it as below. 


2u u*(r) 


3x 


fO 1^ . I . j_- .1^1 


, 1 3T 


or, 


2u 


X 3T 3 i ^ 1 ai . 

u*(r) = — + 7 ^ + 


ar 


1- ifi 

r 30 


Discretizing this equation using finite difference 
approximation for any point (i, j) in the fluid region and at a 
x-station (k-i-l), 




"A"x 


( Ar^)‘ 


k+1 


J I ■ Ft. . - “T. . •! + ■ is t*F. . . 

^ ^i 2^^f ^ 1-1.3 k+1 r? ( ^0)2 

4, A 

- 2T. . + T. . . ] 

1,3 1,3+1 k+1 


For simplification the two dimensional notation for T 

is used, replacing T^^j^k ^i,j,k+1 ^i,3* 

Simplifying the above discretized energy equation for 

fluid and writing it in the standard form of tri-diagonal matrix, 
as below, 


A.T, 


i^i-1,3 


+ B,T. 


1 i,j 


+ C. T 


i*i+1 ,3 



where 


2r7 - Ar- r 
1 f 1 


Bi 


Ar - 2 2 - 4 ^ "I 

= - [4 r^ + 4(-^) + C r^ u^l 


(3.2a) 
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where 





C 


2 

= 2rj^ + Ar^ 
Ar- 2 

= - 2(-^) 

4u(Ar^)^ 




C uf 




Equation (3o2a) is valid for i = 2 to i e= - 1. As 
in the case of solid region, this equation is also modified 
for j =i= 1 and j = and the modification takes place in the 
expression of D . 

The boundary condition (2.4a) is discretized as below. 
For i = 1, i.e. r - 0; ^ 

discretizing it; — Ux_ - o 

or, - T. . + , s= 0 {3.2b) 

* f J J 

Combining equations (3,2a) and (3.2b) the general form of the 
discretized energy equation for fluid is obtained as below, 


A.T. ^ . + B.T. . + C-T. . == D. 

1^1-1, j 1 i,j 1+1,1 1 


(3.2) 


where A^ = 2r^ - 


h 


„ Ar- 2 o 

= - [4r^ + 4(^) + C ] 


C. = 2r^ + ATj 




A^f 2 


th.J-l + + Tp_ _ 

valid for i = 2 to i = - 1 


v- 
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and also for i = 1 , = 


% 


0 

- 1 
1 
0 


Equation (3.2) is modified for j - 1 and j = N 


r 


3.1.3 Disciretization of the Compatibility Conditions at the 
Fluid-Solid Interface 

The compatibility conditions (equation (2.5)) are, 


at r = r. , 


k II 

ar 


= k 


f 3r 


and T ^ = T , 

s 


Energy equation for solid, equation (2.3) 


^ j. X ^ ^ L- £l - 0 

r ar ^2 ^^2 


and energy equation for fluid, equation (2.4) 




aft ^ i il . J_ ifl 

3,2 ^ r 3r g^2 


Refer to Figure 3. 

From Taylor’s series expansion for point (i-1,j), 


A-i,j = 


i,j,f ar 


Solving (■^—7) from the above equation 

a r'^ i , j , f 


(4) 

ar i,j,f 


[T, , , - T, , + Arf(|i), . ] 




At the interface u* = 0, also == • 
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FIG. 3 FINITE DIFFERENCE NODES NEAR THE 
INTERFACE 


i 
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Putting all these in the energy equation for fluid, 
i.e. in equation (2.4), 


0 := 


! \2 L^i-I 1 " ^i i ^ 

(Ar-)^ ^ ^ ^^i,j,f 


[T, 


4 . . 1 

1 i,j,f r^ 30^ i,j,f 


Simplifying it, 


{— ) 
^ar'' 


2r 


1 




2r^ Ar^ + ( Ar^)' 


[T. . - T. , . 

^ 1-1, J 


(Ar^)^ .2 


(ill) 

^ »2l 


2r^ 30^ i,j,f 

(3.3a) 


] 


Expressing T. . . in Taylor’s series as, 

1+1 ,j ,> 


*^i+1,j " ’^i,J 


i ,i , s 


3r"^ i,j,s 


thereby, 

(4) 


ij,s 


”)^ ^^i+Ui ” ^i,j “ 


^T^ j 

i,J,s 


Using this in the energy equation of solid and simplifying it, 

) 

i,3,s 


3‘T 

the (^) term is expressed as below, 


(— ) 

^ 3r''. 


i,J,s 


2r. Ar^ 
1 s 


2r, 

~TTar7P h+i,j u,j 
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-5 3 

2r;j^ 1,0 , s 


(3.3b) 

Ultimately two expressions are obtained for the temperature 
gradient at the interface for fluid and solid regions as, 


(91) 

dr\ . - 

1.0, f 


2r 


1 


(Arf )■ 


Arx(2r, + Ar^) ^^i-1 , j “ — 2 ' ^2 


(=) ] 


V ^ 


2r:j- 0 i,0,f 


and 
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^ » J * ® 


2r 

_ ^ r "T 

Arg(2r^ - '- i+i,j 




(Ar 2^ 

s -■ f g T N 
O ^ O' 

2r^ 30^i,j,s 


The compatibility conditions (equation (2.5)), i.e. Tj^ = 


at r = and 


■ill == k ill 

s 3rl f arU 

s f 


at r = r. 
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( — and (■^— k) will be equal and both of them have 

3ri,j,f 90 i,j,s 

same expression because the temperature of fluid and solid are 

2 

9^T 

the same at the interface and — ^ is expressed as, 

30^ 


(H) 

90 i,j,s 


= (A) 


= (A) 


30 i,j,f 30 i,j,r=r^ 




Using the heat flux equality condition, i.e. 


k (il) 


i,J »s 


-r* 




and expressing (|-p) 3nd (|^) 3s in equations (3.3a) 


, , 3.,j,s i,j,f 

and (3.3b), the following equation is resulted 


k 

s 

IcT * 




1 r r — T . + X — — x) 

f 2r, Ar^ - (Ar^? 2r^ ^i.j 


2r 


1 


2r^ + (Ar^)^ 


l^'^i,j " "^i-l ,j " 


(i!|) ] 




.2 

Using the expression of (-—x) in the above equation and after 

90^ i J 

necessary simplifications the standard discretized equation for 
the interface is obtained as belovy. 



(3.3) 


where, 

B. 


2r — C 

1 “ W 


= - (2r^ R + 2r^ + RR) 




where R = 


and 


RR = 


2r^ Ar^ + (Ar^)^ 

* ‘ ' ■ » ‘I ■'■ ^■■■0.^ ■•mil ■ m m 

2r^ Ar^ - (Ar^)^ 
R ( + ( Ar^)^ 

u7)^ 


Equation (3.3) is valid for all j’s except j == T and j = 

For these two j-stations the equation is modified by replacing 
(j-1) by for j = 1 and { j+1 ) by 1 for j := 

3*1 *4 Simplified Equations for and (x. 0) 


27^ 

/ / u T r dr dp 

From equation (2.6), T^. = , as 

^ ^ 271 

/ / u r dr dp 

o o 

PCp = Constant for the fluid. 

Putting the expression of u(r) = 2u.u*(r), where 

^ 2 

u*(r) “ 1 - (—) , then simplifying the above expression fox 
the following expression is obtained, 

2 2'it 

"^fb ~ 2 u*(r) T r dr dp (3.4) 

nr^j o o 


The heat flux, q^(x, p). 
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qi(x, jZ) = -P 


= k. 


Ar 


r=r. , s 


(3.5) 


3.2 Mumerical Solution of the Discreti zed Equations 

The set of discretized equations to be solved are 
listed below. 

Equation (3.2), j 


°i 


and 



= 2r2 . r. 

Ar^ 



Ar 2 

h 

= _ [4r? H- 

1 



= 2r? + 

Ar^ 


Ar. 

2 

Di 

= - 2(^0 ) 

[ T- • 1 + 

’■ 1,3-1 


T. - C u? rf Tp 

for i — 2 to i = ^rf ” ^ 


^ ° 

- 1 


Gi = 1 


D. - 0 

1 


for i = 1 


Equation (3.3), "^i^i-ljj 


°i 


where \ 


h 


= 2r - C. 

1 


= - [2r^R + 2r^ + RR] 


, i 

K-’k 
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^i 




R 

^3 2r^ Ar^ + ( Ar^)^ 

‘‘f 2 r, Ar^ - ( Ar^)2 


and 

RR 

R( Ar^)^ + ( Ar^ 

( A0)2 


Equati ons 

( 3 . 1 ) and (3.1a) combinedly written as 





where 


= 2r? - r. Ax^ 




0 AT 2 

= - t4ri + 4(^) ] 

for i 


=1 

= 2r? + r. Ar 

1 1 $ 

to i 


h 

Ar 2 

*-^i,j-1 ^iJ + 1 ^ 


and 

\ 

= ar^ 



Bi 


2 

] 



= 0 



D. 

1 




(3.1b) 


^rf ^ 
- 1 


J -3 


where S^ = Arg( 2 r 2 + r 2 valid for i - 


All the above scfuations are to be modified for j = 1 and j s: N^, 
the modifi cations will be only in the expression of . For 
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j = 1 , (j -■ 1 ) is to be replaced by and for j = ( j + 1 ) 

is substituted by 1. 


<r 


f ‘i 
J 


k 

J f 



Equation 

( 3 . 4 ), 


2 

^ / f 

icr ^ o o 

1 

U ' 

Equation 

( 3 . 5 ), 

q ^( x , 0 ) 


J 

- k 

r 

Equation 

( 2 . 7 a ), 

( x ) 

2ti 

dp 

Equation 

( 2 . 8 ), 

f ^( x ) 

- 2ti 

2ti ^ 1 

0 

dp 

Equation 

( 2 . 9 ), 


q - jCx ) 



- ■ffb 


Equation 

( 2 . 10 ), 

Nu 

X 

^ ( 2 rp 





Vf’ 


The set of discretized ec(uations (3.2), (3.3) and (3.1b) for 
fluid, interface and solid are to be solved by a suitable nume- 
rical technique. 

In the present analysis a line-by-line method, sugges- 
ted by Patankar [lol, is used. The line-by-line method helps 
to solve the set of simultaneous linear algebraic equations by 
Tri-Diagonal Matrix Algorithm (TDMA) method. 

The line-by-line method is described in Appendix-B in 

great detail. 
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3.3 Algorithm for the Humerical Solution 

To solve the discretized equation (3.2), (3.3) and 
(3.1b) to get the temperature at a particular x-location, the 
following steps are taken. 

(a) The computations are stated at the location x = 
using the known inlet temperature at the entrance of 
the tube. 

(b) The redirection lines are selected to initiate the 
solution procedure by line-by-line method. 

(c) All the temperatures are assumed initially, both in 
the fluid and solid region. 

(d) Start solving by TDMA with the line = 0. 

The temperature at all the grids of this line are solved. 

(e) Go to the next line along 0-direction (i.e. at 0 - ^0) 
and solve the discretized equations to get temperature 
at all nodes on that line. 

(f) This procedure is followed one after another for all 
the r-direction lines. 

(g) The steps ’d’ to ' f ’ are repeated till the process 
converges, i.e. the temperature at any node does not 
change further with iterative procedure. 

After getting the converged values of temperature, 
they are stored for that particular x-statlon. Then the fluid 
bulk temperature, the average heat flux, the local values of 
. convective heat transfer coefficient and the Nusselt number 
are calculated following the equations (3.4), (3.5) and (2.7a) 


e 
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to (2.10)j The integrations involved are done by Simpson's 
rule of numerical integration. 

The procedure then marches in x-direction and continued 
till the fully developed situation is reached. The fully deve- 
loped situation is indicated when there is no considerable 
change in Nu^ along x. 

The accompanying Flow-chart (Fig. 4) represents the 
algorithm pictorially. The input data to the program and the 
listing of the computer program are given in Appendix-C and 
Appendix-D respectively. 



^ START^ 


^ — ' - - ' 

. 

Initialise all the input data to the program 




Calculate the configuration factor at different 
yJ-location by calling the subroutine 'CONFIG' 



■ ■ ■■ " ■" ' 

k 1 




Guess the temperature of all nodes in (r, 0) 
direction at this x-location to initiate the 
iterative process 

.. . 


t . 


Solve the temperature of all nodes in r-di recti 
at by calling subroutine * TRIDAG’ and in 

this way ail the '^-constant* lines are solved 



STOP 


FIG. 


4 FI 
VU 
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CHAPTER 4 


RESULTS AND DISCUSSIONS 


Figure 5 shows the temperature distribution in the 
tube wall for different conductivities of the tube material at 
the circumferential location, 0=0® (highest point on the 
tube). The temperatures in the tube wall and also the temper- 
ature gradient are higher for low conductivity material than 
for high conductivity material. This is because for low conduc 
tivity material the resistance to heat transfer is high and 
hence higher temperature-gradient. Also, because of less heat 
transfer to the fluid, the inner wall is cooled less by the 
flowing fluid which results in raising the temperature of the 
wall. This suggests that in case of low conductivity tube 
material there is a chance of burn out of the tube at the poin1 
of highest temperature which in turn results in leakage of the 
fluid. 

In Figure 6, the effect of wall thickness on the temp- 
erature distribution in the tube wall is depicted. It is 
observed that with the decrease of tube thickness, the temper- 
ature gradient decreases and the temperature level in the tube 
increases. So, a thin tube with low conductivity is not 
desirable as it may give rise to hot spots occurring due to 
high temperatures on the surface of the tube. The nature of 
the curves satisfy the physics as with the reduction of thick- 
ness, the thermal resistance to heat flow decreases and hence 



TEMPERA 



RADIAL VARIATION OF THE TEMPERATUF 
TUBE VI^ALL FOR VARIOUS CONDUCTIVITIE 
TUBE MATERIAL AT =0* 


(d = 0-1 ,Th = 0-02,W = 10,Pr = 51,Re= 2000 
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FIG. 6 


NODAL POINTS IN THE TUBE WALL 

TEMPERATURE VARIATION IN THE TU! 
FOR DIFFERENT THICKNESS OF THE 
(d = 0.1,Wr1-0,ks =48^Pr=5 
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decreases in temperature gradient. Also most of the radiation 
heat input is used to raise the temperature. 

Figure 7 shows the effect of vertical location (d) of 
the tube with respect to the heated wall on the temperature 
distribution in the tube wall. As expected, the temperature 
level shoots up when the tube is placed closer to the wall as 
it receives more heat from the wall. However, the temperature 
gradient remains independent of the location of the tube (d) 
as the tube thickness and thermal conductivity are same for 
all these cases. 

Figure 8 shows the circumferential inner tube wall 
temperature distribution (fluid-solid interface temperature) 
at a particular x-location for different widths (W) of the 
heated wall when the tube is placed equidistant from the two 
edges of the wall. As one would expect, the curves show symn 
trie nature about = 180o as the radiation heat flux betweer 
Oo and 1800 ig same as that between 180® and 360®. However, 
the curve for bigger width of the wall is placed above that 
for smaller widths of the wall* Also it is more flattened. 
This is quite expected, as the tube receives more radiation 
when the heated wall is more wide* The flattened nature is 
because the circumferential radiation heat flux is more unif 
in nature in case of larger width (W) of the heated wall. 

Figure 9 indicates the plots of mean inner wall 
temperature against the axial distance (x) for different 
Prandtl number fluids, for a constant Reynolds number. It if 
observed that the mean wall temperature is lower for higher 



TEMPERATURE (“K) 


I 

t 



FIS. 7 RADIAL TEMPERATURE DISTRIBUTION IN THE TUBE 
WALL FOR DIFFERENT VERTICAL POSITION OF THE 
TUBE (d ) WITH RESPECT TO THE HEATED WALL 
(Th =0 02 , W = 10 ,ks =48, Pr =51 ,Re = 2000 ) 





i> (DEGREES 3 


^ OF THE FLUiD- SOLID INTERFACE TEM 
HS OF THE HEATED WALL , AT x = 0-1 r 
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AXIAL DISTANCE, X (m) 


FIG.9 AXIAL DISTRIBUTION OF MEAN WALL TEMPERATURE 
DIFFERENT PRANDTL NUMBER (Pr) AT A CONSTAh 
REYNOLDS NUMBER 
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Prandtl number fluid and the temperature gradient near the 
entrance of the tube is more in case of low Prandtl number 
fluid. This is because for low Prandtl number fluid, the 
thermal diffusion is much more as compared with high Prandtl 
number fluids. But, as the axial distance increases i.e. at 
larger x, the temperature gradient is not affected much by 
Prandtl number as the boundary layer thickness grows along 
the x-direction. 

Figure 10 shows the local Nusselt number, varia- 

tion along the axial distance, x for various Prandtl number 
fluids, Reynolds number remaining same* As expected, Nusselt 
number is very high at the entrance of the tube and then it 
decreases monotonically along x until it becomes almost cons- 
tant far downstream. Local Nusselt numbers are less for 
lower Prandtl number fluids. This is because lower Prandtl 
number fluid have higher thermal diffusivity and hence convec- 
tive heat transfer Is less. So, the heat transfer coefficient 
is also less which in turn implies that Nusselt number is less 
Also, it is seen that thermal entry lengths for higher Prandtl 
number fluids are much more as they have lower thermal diffu- 
sivity which means that thermal boundary layer develops slowex 
Figure 1 1 shows how local Nusselt number, Nu^ changes 
along 'x' for various Reynolds number, Prandtl number rema- 
ining same. For higher Reynolds number of the flow, Nusselt 
numbers are larger. This is because, higher Reynolds number 
means higher average velocity of the flow and hence more heat 
is carried away by the fluid and since Nusselt number is the 
indicator of heat transfer, this results in higher Nusselt nu. 




3 20 40 60 


AXIAL OISTANCE,x (m ) 

0 AXIAL DISTRIBUTION OF LOCAL NUSSELT NUMBER FOR VARIOUS 
PRANDTL NUMBER (Pr) AT CONSTANT Re 
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Figure 12 shows plots of fluid bulk temperature 
along the axial distance (x) for various Reynolds numbers, 
Prandtl number remaining same. The plots indicates that 
fluid bulk ■ temperature increases with the axial distance, x. 
For higher Reynolds number of the flow, the fluid bulk tempe: 
atures are less than that for lower Reynolds number. This i* 
because for higher Reynolds number, the average velocity of 
the flow is higher which results in more convective heat trar 
fer giving rise to lower bulk temperature of the fluid at a 
given x-location. 

Figure 13 shows plots of fluid bulk temperature (T^j^ 
along the axial distance (x) for various Prandtl number 
fluids, Reynolds number of flow remaining same. The plots 
indicate that fluid bulk temperature increases with axial 
distance, x. Also, for higher Prandtl number, fluid bulk 
temperatures are lower. This is because for higher Prandtl 
number fluids thermal diffusivity is lower which in turn 
results in higher convective heat transfer and hence lower 

fluid bulk temperature at a given x-location. 

Finally, the accuracy of the present numerical solu- 
tion is checked by comparing the Nusselt number for the full 


Jeveloped condition when constant heat flux condition is 
imposed on the periphery of the tube and conduction in the 
tube wall is neglected, l.e. the tube is considered very ver 
thin. Musselt number at a distance of 2000 meters along the 
tube is found to be 4.5 which compares very well with 4.364, 
the Nusselt number for fully developed constant heat flux 



AXIAL DISTANCE,x(m) 

AXIAL DISTRIBUTION OF FLUID BULK T 
TURE FOR VARIOUS VALUES OF Re OF 





FIG. 13 AXIAL DISTRIBUTION OF THE FLUID BULI 
TURE FOR DIFFERENT VALUES OF PRANI 
( Pr ) 
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condition as obtained by analytical solution. Since, theore- 
tically thermally fully developed condition occurs at infinity, 
it is not numerically possible to go to that extent because 
of limits on computation time. As the local Nusselt number 
for uniform heat flux condition obtained numerically is very 
much closer to the value obtained by analytical solution, the 
present solution for non-uniform heat flux condition can be 
called accurate enough. 
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CHAPTER 5 

CONCLUSIONS 


A simple but computationally efficient numerical 
scheme has been used to solve the problem of conjugate heat 
transfer in a laminar flow of high Prandtl number fluid in a 
circular tube whose periphery is exposed to non-uniform 
radiation heat flux from a large heated wall. The numerical 
results obtained agree very well with the physical expecta- 
tions. Also, the accuracy of the present numerical solution 
has been checked by comparing the Nusselt number obtained 
numerically for thermally fully developed condition with that 
obtained analytically for the idealised case of constant heat 
flux condition on the periphery of a very thin tube. Also, 
this analysis can be modified to take into account variable 
property fluid and the tube wall. The present analysis is 
for infinite length of the tube, but it can also be extended 
to the tube of finite length with necessary change in boundary 

conditions. 


)■ 

tL 
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appendix a 

determination of the configuration factor 

A . 1 Definition 

The configuration factor is also called the geometric 
view factor. For any two given surfaces, the orientation 
between them affects the fraction of the radiation energy 
leaving one surface that strikes the other surface directly. 
Therefore, the orientation of the surfaces plays an important 
role in radiation heat exchange. The physical significance oi 
the view factor between two surfaces is that it represents the 
fraction of radiative energy leaving one surface that strikes 
the other surface directly. 

A. 2 Hottel* s Cross-String Method [ll] 

Consider an enclosure, as shown in Figure A— 1 , consis- 
ting of four surfaces which are very long in the direction 
perpendicular to the plane of the figure. The surfaces can 
be flat, convex or concave* To find the view factor, 
imaginary strings are assumed, shown by dashed line in the 
figure, they are tightly stretched among the four corners A, 

B, C, D of the enclosure. If (i - 1 , 2^3^4,5 , 6 )denote the 
lengths of the strings joining the corners, as illustrated in 
Figure A-1 . According to Hottel, the view factor ^ ^_2 
expressed as, 
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A 


FIG. A-1 DETERMINATION OF CONFIGURATION FACTOR 
BETWEEN THE SURFACES A, AND Aa BY 
HOTTEL’S CROSS-STRING METHOD 


I 



^ crossed length -■ Total uncrossed leriqth 

2 X Length of surface 'V 

“ 51 ; — , Length of surface'!* 

1 ' approximated to 

deriva tion of the Equations for Determining Configuration 
t* 3C to J 

For the present case the configuration factor is 
found between any point of the outer surface of the tube and 
the large heated wall placed over the tube. Each point on the 
upper surface of the tube receives radiation heat flux from 
the entire plate. 

Figure A-2 is considered, the tube is placed at a 
distance d vertically below the large wall. Both the wall 
and the tube are infinitely long in the direction perpendicular 
to the plane of the paper. 

In Figure A-2, A is point of tangent on the outer 
circle from the right end point N, of the wall. Similarly A' 
is the touching point of the tangent drawn on the outer circle 

from the left end point, N' , of the wall. 

Clearly the region below the points A and A' do not 
get any radiation heat flux directly from the hot plate. 

Thereby the configuration factor in this region is obviously 
zero. So, the location of points A and A' is most Important, 
however it is to be noted that when the tube is placed verti- 
cally below any edge then point A or A' will be shifted to 
point C or C respectively. Therefore in the Figure A-2 
angle a cannot be less than P and limiting value of (t is p, 
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however if the top point of the tube, P is placed not directly 
below the plate then a can be less than p. 

To find out P the triangle MCO is considered, 

MO - d + r^ 

OC = r2 “ Outer radius 

,* mg = - OC^, as MCO is a right angle triangle and 

ImO = 900 


or MC ^ /(d + r^)^ -^2 = + 2d r^ 

V^d^ + 2d 

Asj Z.M0C = P , *’• tan P = ^ = 3 


P 


tan 




•f 2d r 
d 


(A.1 ) 


Consider the triangle OAB, 

ZAOB = (1800 -a) 

/.OAB = 900 

9^ = 1800 [/.AOB + /.OAB] = a _ 90® 

In the right angle triangle OAB 

QA _ 9 = sin(a - 90®) = - cosa 

OB I 

OB ^ - OA/coso£ = - r^/cosa 

mb = MO + QB = d + r^- seca 
MN = W - L 
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or, 


or, 


M _ 4. ^ 

- tan = tan (a - 90o) = - cotc 

^'7 - L 

^ + sec a 

('^7 - L)sin ct = „ (d + r 2 )cosa + r^ 


Squaring both sides of the above expression and simplifying it 
the following expression is obtained 
o 

A cos a - B cosot + C = 0 , 

where A == (W - L)^ + (d + r^)^ 


B 

C 

cos a = 


= 2r^(d + r^) 


- (W - L)^ 


B + - 4AC 

2k 


B ± D 
2A 


D = - 4 AC 


sina 


= '/’l - cos^a = ■i^/4A^-(B + D) 


B + D 

when, cosa = — 2 S — 


or, 


= ^V'4a2 - (B - D)^ 

B - D 

when, cosa = — 


tana - 


^4A^ - (B - 'QY 
B - D 


(A.2) 


V4A^ - (B + P}'^ 

B + 9 
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Similarly, to find out 6, the triangle OA'B* is considered. 
In triangle OA’ B’ , 

/.OB’ A’ = ©2 » /.B'OA' = (1800 _ 6) 

= 1800 _ [900 + (1800 _ 6 )] = (6 - 900 ) 

OB' = OA'/sin 0^ = - r^/cosd 

From triangle, MB'N', 

MB' = MO + OB' = d + r^ - sec6 

MN' = L 


tan 0, 


MN' 

MB' 


Td^Tr^-r^secfiJ 


or, - cot 6 - <j + - r 2 sec6 


Solving this equation for '6* 


tan <5 = 


or. 


Aa^ - 

Aa^ - (B^ - )' 


(A.3) 


where, == '^(B^ - ) 


= + (d + r«)' 


B = 2r2(d + T^) 


r2 

I? 2 -■ Lf 


Y 


Also, 


211-6 


(A. 4) 



5d 


From the ecfuations (A.1), (A. 2), (A. 3), (A. 4) the values of p , 
a, b and y are obtained, which play important role for the 
calculation of configuration factor. 

Next step is to express the configuration factor, 

^2W ^ function of '0', Refer to Figure A-3. 

Let, the configuration factor for the elemental arc XY 
on the upper surface is required to be found. The length XY 
is sufficiently small such that in the limiting case this 
element turns to a point and the expression for this 

element becomes the value of at the point P. 

Points X and Y are joined with points M and N by four 
strings, two crossed, two uncrossed, as illustrated in the 
Figure Ar-3. Here, XN and YN* are the uncrossed strings and 
XN' and YN are the crossed strings. 

By Hottel’s Gross-String method, as described earlier, 

at P = 0 

= Lt 
XY-* 0 



From the Figure A-3, XN* - YN’ — XY 

and YN - XN SJ YZ 

where XZ is the perpendicular drawn on YN from X 
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FIG. A -3 
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For 0 > cXj F2 ^ji| - 0, till 0 = y in the other half. 

To find out YZ, the enlarged view of the triangle XY2 is 
considered, as shown in Figure A— 4« 

From this figure, , ^YPO = 90o as OS* is 

normal to the surface at point P. 

Xy is very small, so that YN and PN are almost para- 
llel, and as XZ is perpendicular to YN, so XZ is also perpendi- 
cular to PN, i.e. /.PX^X = 90o in the triangle PX^X and as OS^ 
is perpendicular to XY, /IXPS^ = 90©, 

/.NPX = 90© - 

Zx^XP = 180© - [/NPX +ZPX^X]. = 


Now from the right angle triangle YZX, 


YZ 


^ = sin i.e. YZ = XY sin 

+ ^ 


YY , y 7 ^Y sin Tj 

^2W at P ^ YY^ n 2X7 ^ n 2XY 

Ai U XY -*■ 0 




= -5(1+ sin y^), for 0 ^ 0 < a 


(A.5) 


is found out considering Figure A-4, 


From triangle PS^N, 

From triangle PS2N, 

QM = d + r2, MS2 
00 ^ " ^2 ^2^ 


y^ = 7r - 




tan 



= r2 sin0, MN 
= MN - MS2 = 


+ 11/2 + 0) 
+ 0 ) 


= (W - L), 
(W - L) - 


(A.6) 


sin0 
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FIG. A- 4 
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^^2 ~ OtA - 00^ = (d + r^) - cos0 

(d + r^) - COS0 

•■• = W - L) - r/slni; 



tan”^ 


^ (d + r^) - COS0 ^ 
CW - L/ - sin0 ^ 


(A.7) 


For the left half, i.e. for y < ^ Figures A-3 

and A-5 are considered. 

The uncrossed strings are X‘N' and Y’N, and the 
crossed strings are X'N and Y*N*, all are joined clockwise as 
in the case of right half of the tube. 

As 0 is measured in the clockwise direction, the 
relative positions of points X' and Y' must be noted carefully. 


u T+ (X'N + Y«N» ) > (X>N* + Y'M) 

^2W at P' "x'Y'«*0 


Using X'N - Y'N = X'Y’, and X'N’ - Y'N’ ~ X'Z' 


2W at P 


I — 


Lt 

X'Y' 0 


X'Y' - X'Z' 

2X'Y' 


Figure A-5 is considered to find the expression of 


As ^Y’P'S 3 = 90®. Zy’P’Z; = 90® - 

Again Y'Z' is the normal drawn on X'N’ and for a very small 
element X'Y* , P'N' and X'N' are almost parallel. This gives 


ZP'Z’Y' = 90® 

/.Z'Y'P' = (ZY'P'Zj +Z.P'Z^Y') = ’*^2 

In the right angle triangle Y'X'Z' , 


/X' ^ 90® - ^2 

X'Z' “ X'Y* cos(90® — ^2^ ~ 




X'Y' sin ^2 
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From expression of at P’ , 


2W at P' 


T . X'Y’ - X'Y' sin 'fo 
Lim , - ^ 

X'Y’ - 0 


^21V at P' = = '5 ^2^ 

From the same figure, /OP’O^ - (0 - ir) 


(A. 8) 


M' P'N' := /.OP'N’ - Z0P'02 

= (it - ’i^) - C0 - = 2it - (0 + '^2.^ 

2. ^ 

From triangle N’M’P', M’P'N’ = | *- ©2 


Equating these two, 

2it - (0 + ^2^ = i - ®2 

or, - 0 


(A.9) 


where, ©2 “ ^ ^ 

P'M' = U'O^ + + ^2 ^2 

= d + r 2 - 3:2 COS 0 

M’N' = MN' - MM' = L - r 2 sin{0 ~ %) = L -f r 2 sin0 


^ d 4 r2 >- r2 COS0. 


«2 ^ L 4 To H iT 


(A. 10) 


Therefore all the expressions 
of 0 are listed below. 


for Fow(0) for different values 


4 
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= ■5(1 + sin for 0 ^ 


= 0 


for a < 


= ^(1 - sin ^2^ Y S 


where 


f 


1 “ I - (^1 + ^), 


and Q. - tan 


-1 r ^2^ ~ ^2 1 

L /111 T % 'T" 7 ' Z/i -J 


'1 - "(W - L) - r„ sinjZ) 


= 1.511 +0 0, 


. (d + r«) - r^ cosj2 ^ 
r L V-" r — rTTfli -1 


®2 ^ ^ L + r2 sinJZl 


0 < a 

0 < Y 

0 < 2ir 


and 


o4 


appendix b 

THE LINE-BY-LINE METHOD 


A line-by-line method is a conventional combination 
of TDMA (Tri-Diagonal Matrix Algorithm) method for one dimen- 
sional situations and the Gauss-Seidel iterative method. Acco 
ding to this method a grid line (in a particular direction) is 
chosen and the values of temperatures, T*s (or any dependent 
variable) along the neighbouring lines are assumed or known by 
their latest values, and the T' s along the chosen line are 
solved by TDMA. This procedure is followed for all the lines 
in the same direction and then the technique is repeated till 
the required convergence is attained. Sometimes to make the 
convergence faster the same procedure is repeated for the 
lines in other direction. 

To visualize the line-by-line method more clearly the 
Figure B-1 is considered. The discretization equations for th 
grid points along a particular line (either in x-direction or 
in y-direction) are considered. These equations contain the 
temperatures at the grid points (shown by crosses) along the 
two neighbouring lines. If these temperatures are substituted 
from their latest values, the equations for the grid points 
(shown by dots) along the chosen line would look li^e one 
dimensional equations and could be solved by TDMA. This pro 
cedure is carried out for all the lines in y-directlon and 
may be followed by a similar treatment for x-dlrection. 
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The convergence of the line-by-line method is faster 
because the boundary condition information from the ends of 
the line is transmitted at once to the interior of the calcu- 
lation domain, no matter how many grid points lie along the 
line. The sweep direction (i.e* the sequence in which lines 
are chosen) plays an important role in the convergence of the 
iterative technique. A sweep from upstream to downstream 
would produce much faster convergence than a sweep against 
the stream. 
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APPENDIX C 

INPUT DATA TO THE PROGRAM 


i) Geometry of the tube; 


M 

0.12 


Th = (r^, “ ), Cm) 

0.005 

0.-01 

0.02 

0.03 


ii ) Position of the tube; 


d (m) 


W (m) 


L (m) 


0.07 
0.10 
0.1 3 


0.4 

1,0 

20.0 


w 

2 


, w \ . 

iii ) Conductivity of the tube material, 
k = 35, 48 and 55. 

iv) Known temperatures (®K) 

T = 420.0, Tg = 300.0, = 830.0 

v) Radiation parameters: 

a 0.567 X 10"”^ W/Cm^.K"^), ^ 

vi ) Flow parameter; 

Re = 100, 1000, 1500 and 2000 

where, n (2r^) 

Re is defined as: Re — y 

•-i taken for numerical solution: 

vii) Number of grid points taxen lo 

Nj, = 50, Njf = 30, = '‘0 
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viii ) Properties of fluids^ 


s r 


t 


Pr : 

Kf 

a 

1 

* 

t 

V 


W 

m.K 

m 

s 

1 

1 

1 

1 

% 

t 

2 

m 

s 



' 7 

[ X 10^ 

! 

1 

1 

X 10^ 


1 

t ' 

t 

1 

\ 

L- 


51 

0.260 

0.932 


4.75 

84 

0.132 

0.663 


5.60 

175 

0.135 

0.710 


12.40 

490 

0.138 

0.769 


37.50 
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APPENDIX D 


LISTING OF THE COMPUTER PROGRAM 


NOMENCLATURE X 


A 

~ coefficient 

of 

TDM ! 

NR 

= 


AK 

1! 


P 

1 

NRF 


«rf 

AKF 

k^ 


t 

J 

NPR 

= 

Pr 

AL 

= L 


1 

f 

PI 


IT 

ALPHA 

= <jc 


1 

1 

QB 

= 

h 

ANU 

= Nu 


1 

R1 

tS! 


ANUF 

=s i> 


1 

1 

R2 

ss 


B 

= coefficient 

of 

TDM 

RENOLD 

ss 

Re 

C 

= coefficient 

of 

TDM i 

RI 

ss 


D 

= coefficient 

of 

TDM ! 

SIG 

s: 

<3" 

D4 

== d 


1 

T 

= 

T 

DPMI 

= AtfJ 


1 

1 

TFB 

=s 


DR 

= 


f 

1 

TFP 

= 

■^p 

DRF 

* Ap 

J 


t 

J 

TS 

B2 

Ts 

DX 

- A 'A 


1 

TW 

US 


EPS 

=: e 


1 

1 

TUB 

=z 

t; 

EPSN 

= Accuracy factor for 1 

U 

ss 

u 


convergence check 


F2W 

H 

I 

J 


h 

i 

J 


UAUER = u 

u « final soln* 

vector ofTDMA 

W ^ W 


K 


k 




L 



H. 1 E C H THESIS 

IirLtl C 0 N J U C3 A T E HEAT TRANSFER IN A 
L A M I N A R p I p [: FLOW SUBJECTS Ei lU 
N 0 N ■ U N 1 F fj M K A H I A T I 0 N F R 0 fi A L A F? G E 
H i: A T E B WALL J T H E CASE (J F T H E M A L Y 

u N D ri: g E L 0 p F n r i. o u 

MA[N F'ROGRAM 

)K ^ JK >1< * Jli * 5|< >K ;K 5|i 5K J{< >}; JK 5l( )|< :jc )1< ** jK ?}:*** )R * )K * 5|c * )*c ^ 3K * >K )|i He **** >K >K ** >«c )Fc )i( )»c ** >K 3j£ )>: 
COMMDN/AF<EA J./D4 f R2 y AL y W y r2UI 
C 0 M M UN / AR E A2/ A y B y C y i:i y g 

niMENS TON T(l:10y50) yA(50) yB(50) yC(50) yEKKO) y g ( 50 ) y F2W ( 50 > 
niMENSLON TI'F'CSOyUO) y7 EL ( 50) y S (50 ) y U ( 50 ) yRKSO) yTWALLCSO) 
DIMENSION Q(50) 

OPEN ( UNI T^2 L y LiEO ICE^-' ' DSK ' y FILE=' MASTER * OUT ' ) 

AL -0 * 5 f W^-1 , 0 y Rt==-0. 12? R2===0* 14 y D4=^0 ♦ 1 ? 

AK=4S ♦ 0 ? TS=300 . 0 ? SI6"-0 ♦ 5A7E -07 ? TW=8.J0 . 0 ? EPS=0 ♦ 79 

AKF-0.132yNRF=-=30 

NI-M1 =-NRF-1 

ALPHA=0*(f^A3E-07 

RENOLD -1500 * 0 ?NPR”a4y ANUF~0*05AE-04 
LIAMER-^RENOL D)KO * 5H<ANUr/R1 


TYPE. HcyUAgiL-R 

WR ■[ 1 1" ( 2 1 y 3333 ) NPF( y RKNOLD y ANUF y UAgER 


332 :^ 


F0RMAT(2Xf I6y2XyF64 ly2X»li:8*3y2X?l-“10*4/) 


10 


15' 


riRF'~Rl/FUOAT<NRF 1) 

PI ~3 . 14;l5927 ^ HPHI-- 40 » NR-=‘50 

Dr‘HI---1p;L/FL0AT<NF'l-n ) 
riR =■ C R2“"Ri ) /FLOAT (NK' -NkF ) 

PI::: itcDpr I URFi^DRF 

2t>Kui Aim-itrauR 
rC:-=F.UAK/<P2^AKF) 

RR=^ ( RJfcDR^DR f fiRF^URr ) / ( UPH 1 *DPHT ) 

A1-2.J(JR1*R1 

A3--Al!t:R 

A2::="<A3-i-A1+2t^RR) 
S2^2.3KDR>ttSIG*E:PS/AK 
S2-S2;ic ( 2 ♦ )!cR2)i<R2+R2;icDR ) 

DX=-i *0 

DRDP::sDR/riPHI 

riRDPF="-BRF/r.iPHl 

p:psN“~-o , .1. 
rio to I-'tfNRF 

RI<I)=liRF 

LJ(I): ItO'-KK I )*R1<T )/(R1^Rl> 

DO 15 1 = NKI' i .1. ?NR 
R‘T ( I) “RHliR*C I'-NPR) 

DO 200 J".lyNRl' 

DO 205 J--1 yNPHT 

TFPCly 

CALL CONFIGC lyMPHI) 


205 



40 


t: 


I 6 b‘l 


iOl 


-> 


('ONST-^^f OtDK‘F-')f:BRF''){<UAUL:R/(ALPHA)J(nX) 

xx-o*o 

LUJ B02 K--.1v 50 
DU 10 L--tvNR 
DU 40 J.^lvNPFlI 
T( T V J1 -450* 

11099 ITE.F;: -1,200 
TOI D--=T (Nk'vl ) 

DO 1651. l-=2vNRr-l 

A < I ) = 2 ♦ >KF!;i < i ) mn < i ) -ri < .[ ) )|cdf?f 

C( I) .^2.1{RI(I )>KRI( J)-fRI(I)>KDRF 

B< I ) =:-4, :)cRI ( I ))i«R £ ( D-l »:<cDRDPF'*DRDPF--COKST)kRI (I ) sKRl C I ) 5KU ( I ) 

D < I) •-.= "CONSDKRl ( I ) ( X ) JKU < I ) !j<TFP ( 1 v 1 ) ”2 . ^DRDPFsttDRDPDIc ( T ( .£ v 2 1 -f 

ITdvNPHD) 

A< 1)-0.0 

B<1) =•• L*0?CC1 )-=l *0,0(1) “0,0 
DOlOll-NRF-fX vNR 

A( I > = 2 + ^R;t.( I ijkrkiv-p.'K l)*Dk ^ 

u ( I) -2 , )KRi < I ) mi I < .i: ) +nR:«Ri ( i ) 

B ( I ) --4 * )KR.[ ( X ) 5KRI < [ ) -- 4 , *IlRDP:*DRDP 

Li ( .r ) -- , )KriRDF'JHDRDP>|{ f ( I v NF'HI ) "2 , *riRDP5KDRDP)KT < I ? 2 ) 

A(NR>-4.>;iR2’KR2 

C(NR)-'0*0 

B C NR >-11 (NR)- 4 + ^S2)|!T (NR)' 1 

D(NR)-D<NR)-02*F21J<l)*TU**.(.-S2!Kl.-F2W(n)*TS**4.-3.*S2 


iJKKNRv 1 )>i<>K4* 

A<NRF) AlvE(NRF) A'^ J C< NRF >-“A3 
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r)CNF<:K);:-'RFt>|cf <NI<F» 7)“F\R)KT<NF<KfNPH3 ) 

CAi L TFarsAGd 
PO 102 I-1»NF.- 

102 T(Ul)=rOC.13 

C )({>|o|()Jt)K'*)K>lf^>k)K>k5K5tc;|;))f)KiC>t:3}c )!<)«* J=2yNPFn-l 

DO 103 J=2)-Nr‘FII-l 
DC) 210 T=2rNRr-.1. 

A C T ) •■= 2 * ^<RI ( I ) JifRI ( T > "R J ( T )*DRr 
C(l) 2*)*cR[C [):{<RI<l)+F?I<I)5ftDRF 

B ( 1 > •= --4 ♦ :t;Rl ( I ) >KRI ( 1 ) ~4 ♦»tDRriPF>JfDRriPF-CONST)f(Fa ( I ) )KRI ( I ) *U ( I ? 

210 D( ,T )=--C0Nf;T3KR] ( D^RT ( I ) *U ( I >>}fTrp ( I y J )“2t)kriRDPF)KDRDPF)tc ( T( X J+1 )4 

3 T < I y J- 1>) 

A(3 )--=o.o;r;<i)^^i*o 
XH 1) “ U 0 5 D ( t ) ==0 + 0 

DO 3 04 T==NRFFlfNF< 

A ( X ) -^2 * >F.RI ( I ) m<L ( I) “RI < I ) tl\R 
C< X )•••■?* >KR] (] ):<^RI(I)+RX<X)3KIiR 
P ( :i ) ^ -4 . )KRI C X ) 5KRT < I ) ~4 . JKliRriP^FdiRDP 
104 D(X )== -2 .)|cDRDP*URIiF->!«X CLi- J-1) -2**DRDP>kIiRDP*T<lTJ+l) 

ACNR)'-=-4*JKR2'-itF^'2 
C(NR)'-0*0 

B ( NR ) ==-"H ( NR ) ' 4 * *S2)K r ( NR y J) *#3 * 

B f NR)~D(NR)-B25i<F2U< J)*TUI:F:*4»-"S2#( 1 t-F2W( J) )>l(TS^5ft4 ♦ -S + JFtSJ 
l>f{ T (NRy J)*5K4» 

A (NRF)™A3 ?B(NRr- ) =A2 ? C (NRF) "A3 
D ( NRF ) ” -RFa ( T t NRF y J 1 3 ^ +T < NRF r J--1 > ) 

CAl-l. TRIDAGCI yNR) 


rO 105 I lyNR 
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105 T(.r? 

10.3 continue: 

DO 215 1^2yNRF~l 
A ( 1 ) -=2 ♦ SXK H 1.) :icFa ( 1 ) •• RI ( I ) ^diRF 
C < I ) , )({R ;C ( T ) T < T ) {-R I< T ) JkBRF 

B ( 1 > .= - 4 * )|!RI ( T ) *R r < [ ) “4 4 >l<i:iRliPF*DRriPF-CnNST*RI < I) 5}cRI ( T ) )KU < I ) 

215 D ( T ) = -CONS T *R [ t: T > jKRT ( 1 i »:ti < f ) )KTFP ( I y NPHI ) ~ 2 4 >t:DRDPF*riRDPF!t: ( T ( I r 1 ) 

l+T(iyNPHI-1 ) ) 

A f 1 ) =0 4 0 ; B ( t ) -= “ 1 4 0 y C < 1) - 1 4 0 f IK 1 ) =0 4 0 
DC) 106 J:^NRF+lyNF< 

A ( 1) =2 4 >KRI ( I ) ^RI ( I ) --k’T < T ) *DR 
C(I)=--2 4>|cRX< I )5l«R,r<I) IRK [)>tcDR 
B ( J ) --4 4 *R :i. ( T ) )KRI ( X ) "4 4 )KDRDP>lcDRDP 
1 06 D ( :i ) -^-2 4 )fcDRDP>l(DRI,iP)KT ( I y NPHT- 1 ) -2 4 ^jDRDP^DRUP*! ( I y 1 ) 

A(NR)=4 4!KR2*R2 
C(NR)--==0 4 0 

B < NR ) =B < NR ) “4 4 >lcS2)KT ( NR y NPHI ) **3 ♦ 

D ( NR ) =D ( NR ) •S2)KF2W ( NPHI ) :F{TU*3K4 ♦ --82* ( 1 4 “F2W ( NPHI > ) ^TS>K3t;4 4 -S25K3 4 * 
lT(NRyNPHI))K)K4 4 
A ( NRF ) ~A 1 y B ( NRF ) ==A2 ? C < NRF) =-A3 
D < NRF ) ~=-RR*:T ( NRF y 1 ) -RR*T ( NRF r NPHI -1 > 

CALL TRIDAG(lyNR) 
rif3 107 1 = 1 y NR 

107 T< I yNPHD^UC 1 ) 

IF f AEiSCTQLD-KNRyl > ) ♦1-1= ♦EPSN)G0 TO 110 

99 CONTINUE 

WRITEC'^lyllS) 


jp- 


US F"0RMAT(2Xr "NOr CONVI-RflUi V/ ) 

UFaTE<23 U77> 

777 F0RMATC2XUINTFR FACE TEMP* / v NPHJ V > 

WRUE(21U79) (KNRFyJ) j J“1 i-NPHI) 

779 F-0RMriT<8(r7l2*3) ) 

WRl rE<21 

COS F ORMAT (2XU rCHF-. TOR TIIF. WHOLE REGION WITH J=t . NPHT .SS ! -'=1 Ur/) 

lHaiH ( 7U*) ( <T(T !» JU l=--t i‘NF-H[ »5) y 1"= I fNR) 

DO 009 L- I yNRI 
DO 889 .P^UNFTIF 
809 TF“r(I y J) =T( I y J) 

C !l F A I T K A N 0 r L U 0 A L 0 U I A T 1 0 N 0 

DO 11 I = UNRI- 
DO 12 J^^l yNPHi; 

12 TT1 ( J)-T(Iy J> 

TTKNFTII + l) = TT.t (1) 

CALL SrMP0N<0yTPlyNPHi:yTT1 yXINTGI.) 

S<I)=XINT61UR1 < [UTK F) 

11 CONTINUE 

171 CALL SlMP3N<0yRi yNF'MI yO*XTNTOL ) 

TF"B=2*)KXINTGL /<FM:f:R1 *R1 ) 

DO 13 J -lyNFT~n 
TWALLCJ )=•■!< NRP»‘J' 

1,:^ Q < J ) = ( T ( NRF-t U J ) -T ( NRF y J ) > 

TUALLINPHIt 1 )~T WALI < 1 ^ 

Q i NPH ! -f 1 ) -0 ( :l- J 

CALL SIMI-'ONCOy IPTyNPIIJ yTWALLyXINTGL) 


TUB XINTGL/TFI 



CALL SXMPSN(0yTri?NPM [yQ.X]N7i;L ) 

an =x.i:ntgi /‘tpi 

H^AKJKQB/CTWn TFEi) 

ANU=2**H>KR.l/AKF 
TYPE :<<yANU 
XX--=XX+DX 

WRITEC21 yl4)XXy TUB y QB y TFB y H y ANIJ 

FORMAT t 2XyFJ3t2y4XyF84 3y4X»F?*2y4Xy rR,3,4XyF7*3y4XvF£J*3) 

COMTINUF 

STOP 

END 

S U B R 0 U T .C N E J T R 1 Li A (3 

(‘THIS SUBROUTINE SOL.UES SIMULTANEOUS LINEAR NON-HOMOOENIOUG EONS* 

BY TRI- DIAGONAL MATRIX SOLUTION) 

^ 5K >|i :t: )K * )K >1< )K 5K * r3< )K ;K >K :j< 5K * 3|< )K >K 3|c >K * >|J 3tc # >ic )i< He ^ Jk * >K >l< )l< * )if >t« Jf: )f: * )K $ )tc it: )K >K 

SUBROUTINE IRTIiAG< II- vL) 

C0MMGN/AREA2/A y B y C y D y U 
DIMENSION PFTA(51 )? GAMMA ( 51 ) 

DIMENSION A (SO) yp (SO) y(XSO) *0(50) yU(50) 

COMPUTE rNTFRMEU riATE ARRAYS BETA AND GAMMA 

BETA( I1')"~B< 1F:i 

GAMMA( IF) ='D( IF) /BI.T A ( IT ) 

IFPl^im 
DO 1 T==-]rT‘1yl. 

BETA( I)“B(I) -A( I )H(C(I~Sl ) /BETA ( 1 --1, ) 

GAMMA(I) end) A( 1 )*GAMHA( I-l ) > -''BLTA( I ) 


7 / 


C COMPUTF FINAL, SULUTIDN yECTOk V 

V(1.„) =--GAHHA(I. ) 

I.AST==L. -IF 

DC) 2 yl.AST 

[-■=L-K 

2 V(I)--=fjAMMA< J ) J tl>/DETA<I) 

RETURN 

END 

L s u E R G u T I N r-: : C 0 N F I 0 

C SUERGCniNE TU CAL.CULAIE CONFIGURATION FACTOR 

SUBROUTINE CONFIG ( IS » NPHI ) 

COMMON/AREA t /lU y R2 v AL. y U y F2W 
DIMENSION r2Ut<U0) 

PT-=-3* 1410927 
rF‘I=^-2.)f:F'I 

D1~SQRT ( D4>lcD4-t-2* *R2)KD4 ) 

BETA==ATAN<Iil/R2) 

B1~24 *R23KLD4+R2) 

Cl-=--R2*li;2- (W-ADXfClJ -AL) 

A 1 == ( D4 +R2 ) Jf: ( D4 f R2 ) t ( U “Al.. ) JK ( W-AL ) 

ni = S0RT(Bl>KBl""4*>l«Al>FCl ) 

D 2 == 4 . A 1 * A 1 - ( D 1 + E 1 > >1^ t B 1 + B 1 ) 



D2:rSQKTCD2) 

ALPHAl “ArAN(D2/LB14Dl ) ) 
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ri2=4*)fcAlJ|5A.L~(Dl- 

D2=SaRT<rC) 

ALPMA2»A rAN< rnv < B] -DX ) ) 

XF ( ALP1-IA2 * LT ,0*0) ALPX^A2•-ALPHA2+P.T 
.IF- ( AiJ“-'MA2 * I..7 . BFTA ) GO T05 
ALr-'HA"ALPHA2 
XiO 7 0 6 

G ALPHA -ALPHA! 

* C;i=F^2)kR2 -Al.ilcAL 

A 1 =. AL * A L F < rf4 fR2 ) * ( n4 + fC2 ) 

li 1 s-SQRT ( B 1 iKB 1 -4 ♦ ij: A I. *C 1 ) 

ri2 -== 4 , Jic A 1 1C A X - ( D 1 4 B 1 ) < B t X B 1 ) 

rt2'=SQRr<D2) 

/ 

ALPl lAl -=ATAN ( Ii2/ ( Bi FBI ) > 
ri2“4*lcAl3f:Al"-(r(l-Bl)>f:(rH-PX ) 

D2™SnRT(D2) 

ALPHA2-A FAN ^ B2/ ( B1 --HX ) > 

IF' ( A LPHA2 ♦ L 7 ♦ 0 ♦ 0 ) ALPHA2 -ALPHA2-f PI 
.[F( ALPHA2.L7 ♦BETA) 00 TO 7 
riF:(„rA-=ALPHA2 
I'A TO 16 

7 BELTA-'-ALPHAl 

1 6 B P I- 1 1 === 7 F- X /FL 0 AT ( N P H I ) 

OAMA'-TPI-riE'LTA 

PHI ~0 ♦ 0 

K-^^1 

10 X=^D44R2 R21CC08XPH.T) >/(W--AL"R2!tfSIN(PHl> ) 






THCTA=--ATAN(X) 





SYl==-0 * 5Wl~ TMETA-PHI 
FlW=^O.fj:«?U.+GIN(GY:l )) 

r- 2 W(K):--l-lU 

UFiITII(21.>t^)ri!IyFlW 
PHI "PH I -f DPI U 
K“K-fl 

IF- (PHI *LE.ALPFm)GD TO 10 
U; 1 lU-'O.O 

wf<:jtE( 23 y Y)r'Fn yriw 
r2U)(K)-="n w 
K=--K-F'l 

p-HI-^PHI-fDPHl 
II--<PHIyLL- ♦GAMA)GO I 0 Kv 
20 X==(li4-I-R2“R23|CCDS<PHE) )/<AU4-R?3|cSTN(PH1 )) 

THETA-AIANCX) 

SYl:-^1.5)|tPI-f IHETA-PHt 
FlW=0*!:i)K<l. -S;[N(SY1.)) 

WRnE(21 y*)PHIfFlW 
F2W(K)=-~F1W 

PHI=PHI-frtPHI 
ir(PHI.LE*TPI)GU 10 20 
RFl URN 
END 

c ******************««*****«*w*»***************************^ 

c S U B R 0 u r J N F : S 1 M P s H 
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I 

THIS <-7UBPR0GRAH INTI-ret^ATCrS ANY tUOL-N FUNGI TON OVER A FINITE 1 

i 

LJMIT USING SJMI-'SnN^S J/S-RD RUl C I 

i 

■i 

SUBROUTINE SIMI 'SN C XMIN y XMAX y Ny DUMMYF ^ XINT (3L ) 
rUMENSiON nUMMYF^^jO) 

(XMAX“XMIN)/FLOAT<N> 

SUM= 0 + 0 

no 4 r-= 2 rN 

IF (Morui y2) )2y2y;5 
SUM^^-SUM f4 + )KDUMMYr < 1 ) 

GO TO 4 

SUM-SUM+2 + )KriUHh1YF < T ) 

CONTINUE 

XIN T GL'=--H/3 + )i; < DUMMYF" < 1 ) +SIJM-FLIUHMYF ( N+1 ) ) 

RETURN 

END 





